Short-Range Order in a Flat Two-Dimensional Fermi Surface 
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We present the two-loop renormalization group (RG) calculations of all the susceptibilities asso- 
ciated with the two-dimensional flat Fermi surface with rounded corners (FS). Our approach follows 
our fermionic field theory RG method presented in detail earlier on. In one loop order our calculation 
reproduce the results obtained previously by other RG schemes. All susceptibilities diverge at some 
energy scale and the antiferromagnetic SDW correlations produce indeed the dominant instability 
in the physical system. In contrast, in two-loop order, for a given initial set of values of coupling 
constant regime only one of the susceptibilities at a time seems to diverge. 
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I. INTRODUCTION 

In the last decades there have been important develop- 
ments in both experiment and theory which resulted in a 
better understanding of the high temperature supercon- 
ductors. One typical characteristic of those compounds 
is their extreme sensitivity with doping. At half-filling 
the cuprates are Mott insulators. However, at very low 
hole doping, as soon as the antiferromagnetic ordering is 
destroyed those materials are turned into the pseudogap 
phasei^. 

From the start Anderson^ argued that the Hubbard 
model (HM), which considers only short-range interac- 
tions, describes qualitatively well the electronic prop- 
erties of the high-Tc superconductors. In the one- 
dimensional case (ID) this model captures the main fea- 
tures displaying the existence of the Luttinger liquid for 
U > 0, spin-charge separation, as well as a Mott insu- 
lating regime. There is no exact solution for the HM in 
D > 1. Moreover successful methods in D = 1 such as 
bosonization& and the Bethe ansatz£ are either inappli- 
cable or simply much too hard to implement in higher 
dimensions. One approach which is equally successful in 
both one and higher dimensions is the renormalization 
group (RG) method. Different RG schemes are already 
available to describe strongly interacting fermions in the 
presence of a Fermi surface (FS) 8 i 9 i 10 ' n i 12 . 

In this work we continue to explore the field theoretical 
RG method. We calculate the renormalized form factors 
of the main physical instabilities of the two-dimensional 
HM up to two-loop order. We compare our results with 
other RG approaches in 2D2*iS. Those results are sum- 
marized in the phase diagram which displays all leading 
and subleading instabilities in that dimension in coupling 
constant space. To reinforce our case we test our method 
in well known grounds and show it reproduces the analo- 
gous phase diagram in the one-dimensional cassia. In ID 
the coupling functions approach fixed point values pro- 
ducing susceptibilities which exhibit power-law behav- 
ior. This reflects the absence of spontaneous symmetry 
breaking and long-range order. We must keep this in 



mind when analyzing the possible existing states in the 
phase diagrams of such low dimensional systems. Here in 
specifying such phase diagrams we find convenient to de- 
fine appropriately symmetrized response functions with 
respect the spin projections of the corresponding particle- 
particle and particle-hole operators. In this way, we cal- 
culate the corresponding charge density wave (CDW), 
the spin density wave (SDW), the singlet superconduc- 
tivity (SSC), and the triplet superconductivity (TSC) 

(2 1) 

one-particle irreducible functions T R ' s which in turn 
define the respective susceptibilities. 

In a ID metal the FS reduces to two Fermi points 
(±/cf)- In 2D, for non-interacting electrons, at half-filling 
the FS is a perfect square. If we go slightly away from 
half-filling, doping the system with holes, all corners get 
rounded and the area inside the FS is reduced from its 
square value. As a result in such conditions we may ne- 
glect the van Hove singularities as well as the Umklapp 
effects. In this 2D scenario a new symmetry arises with 
respect to the sign of momentum component along the 
FS. This will reflect itself in the RG equations for the 
renormalized form factors which we redefine as the charge 
density wave of s and d-types (CDW±), the spin den- 
sity wave of s and d-types (SDW±), the singlet super- 
conductivity of s and d-types (SSC±) and, finally, the 
triplet superconductivity of s and d-types (TSC±). All 
those symmetries are derived from the particle-particle 
and particle-hole operators which have a transfer momen- 
tum q. Our results for the uniform spin and charge sus- 
ceptibilities were presented in a separated publication 14 . 

We showed in another worki^ that for such flat FS as 
soon as we add interactions to the 2D electron gas the 
flow to strong coupling destroys the Landau Fermi liq- 
uid regime. We demonstrated how the self-energy effects 
lead to the nullification of the quasiparticle weight di- 
rectly affecting the RG equations for the renormalized 
coupling functions. In this work we discuss how those 
results also influence directly the response functions of 
the physical system. In particular we reveal that for the 
lightly doped 2D HM, all susceptibilities flow to fixed 
values. We argue that such a behavior suggests the exis- 
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FIG. 1: The 2D flat Fermi surface with rounded corners. We 
divide it into four regions: two of the solid line type and two 
of the dashed line type. The perpendicular regions do not 
mix in our scheme. 



velocity and the energy e is measured with respect the 
chemical potential [i. Notice that this dispersion rela- 
tion depends only on the momenta perpendicular to the 
Fermi surface, where the label a = ± refers to the flat 
sectors at p± = ±kp respectively. In doing so we are 
only considering the flat parts whose nesting vectors are 
Q* = (0, ±2kp). We consider a fixed momentum cutoff 
A which results in the interval kp — A \p±_\ kp + A 
for the perpendicular component in the vicinity of the 
FS. We neglect how the interactions renormalize the FS 
itself since it would be too complex to do otherwise at 
this stage. As a result, the momentum kp will not be 
renormalizcd in our approach and will resume its nonin- 
teracting value. For the same reason, we will also neglect 
the Fermi velocity vp momentum dependence along the 
FS. 

Since we showed in detail all the necessary steps to do 
the renormalization in such a model up to two-loop order, 
from the start, we rewrite the renormalized Lagrangian 
L associated with the 2D flat Fermi surface entirely in 
terms of renormalized fields and couplings 



tence of short-range correlations typical of a spin liquid 
phase. This is in agreement with our already mentioned 
uniform susceptibilities results^*. 

The presentation of our results begins with the deriva- 
tion of the corresponding RG equations for the vertex 
functions and the associated susceptibilities within the 
field theory method. We discuss initially one-loop or- 
der results and show that they reproduce both the par- 
quet and other fermionic RG schemes. Following this 
we present our new two-loop results. We show that the 
SDW and the d-wave superconductivity are indeed the 
leading instabilities in the physical system when the par- 
allel component of the transfer momentum is equal to 
zero. Finally we conclude our work discussing the physi- 
cal meaning of our results. 



II. THE MODEL 

Firstly we consider a 2D lightly doped FS with no van 
Hove singularity as shown in Fig^ For convenience and 
to keep a closer contact with well-known works in one- 
dimensional physics^, we divide the FS into four regions. 
We restrict the momenta at the FS to the flat parts only. 
The contributions which originate in scattering processes 
associated with FS patches which are perpendicular to 
each other are irrelevant in the RG sense. For this reason 
we can select, for example, the patches centered at (0, kp) 
and (0, —kp) and simply neglect the other two. Following 
this we will restrict ourselves to the one-electron states 
labeled by the momenta p± = k y and p\\ — k x . Accord- 
ingly the momenta parallel to the FS is restricted to the 
interval —A ^ pn ^ A, with A being essentially the size 
of the flat patches. In the same way the energy disper- 
sion of the single-particle states is one-dimensional and 
given by e a (p) = Vp (\p±\ — kp) where vp is the Fermi 



L= Z ^R{a)o (P) \- id t - V F (Kl ~ M] ^fl(a)<T (p) 

p,cj,a— ± 



Y[ Z (Pi\\) [g2B$a5$/3~f - giBSorySps] 

p,q,k ck,/3,<5,7 . 

x ^r(+)s (p + q - k ) ^r(_) 7 ( k ) i>R{-)fs (q) i>R(+)<x (p) , 

(2.1) 

following the "g-ology" notation. Here the ^(i) and 
4>r(±) are, respectively, the creation and annihilation op- 
erators for particles located at the ± patches. The cou- 
plings giB and g2B stand for bare backscattering and 
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FIG. 2: The interaction processes in the model and the cor- 
responding Feynman rules for the vertices. The gm and g2R 
couplings stand for the renormalized backscattering and for- 
ward scattering respectively. 
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forward scattering couplings. From now on we will con- 
sider the thermodynamic limit (V — > oo) with all mo- 
menta summations becoming integrals such as J2 P ~^ 
V J d 2 p/(27r) 2 , with h = 1. Finally, the bare coupling 
functions are related to their renormalized associates by 



9iB 



(g iR + Ag iR ) 



(2.2) 



The diagrammatic representations of the correspond- 
ing renormalized forward and backscattering interactions 
are shown schematically in Fig. 2. Here, we do not con- 
sider the Umklapp processes due to the fact that there is 
no superposition between our FS and the correspondingly 
Umklapp surface. In all Feynman diagrams, the non- 
interacting single-particle propagators G^K and G^\ are 
represented by a solid and a dashed line respectively, 
following their association with the corresponding FS 
patches. They are given by 



<?Mp)) 



0(-e a (p)) 



G (0) (p) = 

(a) po ~ e a (p) + iS ' p - e a (p) - id 



(2.3) 



Before we deal directly with the one-particle irre- 

(2 1) 

ducible function T R ' s. In the next section we sum- 
marize the main results of the field theory for the quasi- 
particle weight and the coupling functions as well as the 
respective RG equations. 




FIG. 4: The diagrams for the renormalized four-point vertex 
in the backscattering channel. 



the corresponding counterterm diagram associated with 
the multiplicative fermion field factor Z. The use of an 
appropriate RG prescription for the one-particles Green's 
function at p± = ±kp together with those contributions 
leads to the determination of Z. It then follows that the 
quasiparticle weight satisfies the RG equation 



III. RENORMALIZED COUPLING FUNCTIONS 
AND SELF-ENERGY UP TO TWO LOOPS 

A. Self-energy and quasiparticle weight up to 
two-loops 

In two-loop order another quantity plays an impor- 
tant role in the RG equations for the renormalized cou- 
plings: the quasiparticle weight Z. To calculate Z we 
must determine the renormalized self-energy Figure 
displays the four important Tir diagrams for the de- 
termination of Z up to two-loop order. Those diagrams 
produce logarithmic singularities multiplied by the factor 
(po — vf{\pa_ \ — hp)) which can only be canceled out by 




+ x + ••• 



" a! =7^(ph;^), (3.1) 

where 7 is the anomalous dimension. The full expres- 
sion for 7 is given in our Appendix A. The numerical 
estimation of Z follows our previous work. As we ob- 
serve later the suppression of the quasiparticle weight as 
we strengthen the interactions changes dramatically the 
one-loop scenario. 

B. Two- loop RG equations for the renormalized 
couplings 

In Eq. I|2.1|) we wrote the renormalized Lagrangian 
which will automatically generate renormalized physical 
quantities in perturbation theory at any loop order. As 
a result, applying suitable Feynman rules one can arrive 
at the diagrams shown in the Fig. 0] and Fig. The 
last diagrams in both sets are the counterterms which 
render the theory finite. In this way, we identify the 

renormalized one-particle irreducible r^j' {% = 1,2) such 
that, at the FS, the corresponding renormalized coupling 
functions gm, (pi||,j>2||}P3||;w) are given by 



FIG. 3: The diagrams for the self-energy up to two-loops. 



H {Pi,P2,P3)\fs = -igin(Pi,P2,P3;v) , (3.2) 
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FIG. 5: The diagrams for the renormalized four-point vertex 
in the forward scattering channel. 



Taking into account the associated Z factors for the 
external momenta and the RG conditions dgiB/duj = 
for the bare coupling functions the RG equations for the 
coupling functions follow immediately 



dgw. (pi||,P2||,P3||) 



duj 



zZ^I foil) 9iR{Pi\\,P2\\,Ps\\) 
dAg iR (pi\\,P2\\,P3\\) 



duj 



-. (3.3) 



where i = 1,2. We make use of the numerical estimates 
of those RG equations in the calculations presented in 
this work. Having said that we are now set to implement 
the RG strategy to calculate the renormalized form fac- 
tors which are indeed the linear response functions of the 



system for the fiat FS in 2D. 



IV. RESPONSE FUNCTIONS 

Following the RG strategy to study the pairing and 
density wave instabilities of the system let us add ini- 
tially to the renormalized Lagrangian of the system two 
fictitious infinitesimal external fields hsc (for the pair- 
ing term) and how {tor the density wave) which act es- 
sentially as source fields for the generation of particle- 
particle and particle-hole pairs. That is, we include in 
our renormalized Lagrangian the contributions 



L ext = I [^ 1/2 (fc||)^ 1/2 (<Z|| - h)hf c (q) T** (k, q) 



k,q 



x (k) v>V)/3 (q - k ) + zl/2 (h )^ 1/2 ( fc n - «u ) 

X hf w (q) T*« (k, q) ^ R(+)a (k) V> fl (_ )/3 (k - q) + H.c] 

(4.1) 

where 7? Ba/3 (k, q) (i = SC, DW) is a bare form fac- 
tor which is determined in such a way to incorporate 
the symmetry it is associated with. Notice that in this 
model the external fields are made spin dependent since 
this is important for the definition of some symmetries 
which we will refer next. As we have done in section II 
we will consider again the thermodynamic limit with all 
momenta summations becoming integrals. By means of 
the added Lagrangian L ext we are now able to generate 
the one-particle irreducible functions associated with the 
composite pairing and the composite particle-hole opera- 
tors. Since we are interested in the response functions for 
the density wave and superconductor channels we need 
to define the associated three points generalized Green's 
functions namely G^^ap an d ■ ^ n doing this we 

get back the corresponding Tp^ap an d by cut- 

ting out the external legs of the corresponding Gf^'^'s. 
The Gf^'^'s are given by 



G S (P. <?) = i -^fs — (^ (+)Q (p)M-)/3(P - ^) ex P 



L, 



qaPa 



,jS , 7 <5 
n DW SC 



75 



(4.2) 



and 



G 



H(2,l) 
SCaf3 



(P, «) = ^ fe £ (q) ( ^ (+)a (p)V^ ( _)^(q - P)exp 



qoPo 



, 7 <5 , 7 5 
n DW'"'SC 



(4.3) 



where < ... > stands for f VipVip^ expiS{ip,ip ji ]( ... ) with S being the classical action associated with the 
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renormalized Lagrangian given by Eq. I|2.1|) . 

Following conventional Feynman rules we display the 
diagrams in Fig. for the density wave and supercon- 
ductor channels up to one-loop order respectively. 

(2 1) 

Notice that the T R ' counterterm diagrams cancel ex- 
actly all the corresponding higher order diagrams which 
scale as (rn(0/o;)) n with n > 2. Moreover, there are no 

(2 1) 

other linearly log divergent diagrams in T y R ' '. Conse- 
quently, in our scheme the resulting higher loop effects 

(2 1) 

in higher order contributions for T R are produced by 
the Z factors in our RG equations. Hence, if we con- 
sider the external Lagrangian i|4.1[l we can rewrite the 
bare form factors such that 



tE$ (p. q) = z-"'' z {p\\ , ^)z-'i\v\\ - % , w) 



- Z-V2 



1/2/ 



x^(p||,?||,w)T^(p,q) (4.4a) 



Bct/3 



SC 



(p, (i ) = Z-^ 2 (p h u J )Z- 1 / 2 (q ll -p h u J ) 



Zsc(Pb<lb u ) T sc (P»q) 



(4.4b) 



or, equivalcntly 



^ ( P) q) = Z-^frii.oOZ-Vafrn - g||lW ) 
x r^(p,q) + Ar^(p,q)" 

r/ c ^ ( P ,q) = z- 1 /^,,,^- 1 ^^ 



(4.5a) 



r^(p,q) + AT s ^(p,q) (4.5b) 



The counterterm ( AT^f^ sc ^ ) guarantees the cancela- 
tion of the divergent vertex functions diagrams in one- 
loop order. This renders the theory finite at the FS. 




FIG. 6: The Feynman diagrams up to one-loop order for the 
renormalized T D Sf in the density wave channel and the renor- 
malized T^q 13 in the superconductor channel. 



Notice that the renormalized form factors depend on the 
momenta p and q. Now, we are ready to set up the pre- 
scriptions for the renormalized one-particle irreducible 

(2 1). 

Green's functions r)^' Q ^'s, in terms of experimentally ob- 
servable physical quantities 



^DWap&bPo = ^,p^ = k F ;q\\,q± = 2k F ) 



^SCaP&bPo =u,p ± = k F ;q\\,q± = 0) 



iT Dw(Pb<l\\) ( 4 - 6a ) 
0) = 

* T Sc a/3 (P||>9||)- 

(4.6b) 



As we already mentioned all renormalized quantities 
are RG energy scale dependent. However, to avoid over- 
loading notations we omit this dependence from now on 
except when it is strictly necessary to do so. In this way, 
considering the diagrams shown in FiglBJand making use 
of the RG conditions stated in Eqs. I|4.6a[) and (|4.6b|) we 
arrive at 



L ^ 1 DW 



4tt 2 v f 1 



9iii(k\\,P\\ - q\\,p\\) 

<r=U 



x TDw{kp1\\)~ 92B.(k\\,P\\ - q\\,P\\) 



(4.7a) 



AT SC P = ~ dk \\ [92R{k\\ , <?|| - fey , q\\ - p\\ ) 

e(*||.9|| - fc lhP||)] r sc Q/3 ( fc lh9||) ln (^) 

(4.7b) 



9ir( 



Consistently with the renormalized density wave and 
pairing form factors and we can now do the 

symmetrization with respect the spin components to de- 
fine 



r cDW^Pb1\\) 



r^(p|,,g||)+T^(p|, J? || 



R 



SDW 
-R 



(p\\,q\\) = r^( Pll ,q ll )-r^(p hqil ) 



Rll, 



T« c (p„ , g„) = (p„ , 9 || ) - T$ 1 (p\\ , n ) 
Tfsc(P\\ , <?||) = T S R c l (P\\ , <?||) + T^(p h g,|) 



Rll, 



(4.8a) 
(4.8b) 
(4.8c) 
(4.8d) 



where (CDW) relates to charge density wave, (SDW) 
to spin density wave, (SSC) to singlet superconductiv- 
ity and (TSC) to triplet superconductivity respectively. 
Now, we can differentiate the equations (|4.5a|) and l|4.5b[) 
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with respect to u to get 



X [ T Dw(P\\ >Q\\) + AT DW Oil > ?|| 

d 



Ra/3 



ui—Tc 
duj 



Ba/3 
SC 



duj 



(jsc P (PU.?||) + AT sc a/3 (P|h9||)) 



(4.9a) 

-^(Pll^)) (V^-^o,)) 



(4.9b) 



Taking into account the fact that the bare quantities do 
not know anything about the RG scale we can arrive at 



5«(Pl|hP2||,P3|hP4||) = 0ifl(-2>l||,-P2||,-P3||,-P4||) 
5«(Pl|hP2||-P3|hP4||) = 9iR{P2\\,Pl\\,Vi\\,PZ\\) 
9iij(Pl||,P2||,P3||,P4||) = 0ifl(P4||)P3||,2>2||iJ>l||) (4-12) 

These conditions need to be satisfied by the RG equa- 
tions in order for them to produce the correct numerical 
results when we approach the FS. These symmetries are 
satisfied by the RG equations (|4.11a(l and (|4.11b(l which 
are symmetrical with respect to the sign reversal of p,, 
for a fixed qu . Following this we can therefore define two 
irreducible representations of this symmetry which never 
mix with each other 



du> 
1 



Tjowtebn) [7(P||.w)+7(P||-g||.w)] (4.10a) 



+ 2 T sc(P||,9y) [7(P||,w)+7(9|| -Pb")] (4.10b) 

where 7 is the anomalous dimension. Since we have al- 
ready defined the symmetrized renormalized form factors 
with respect to spin projection, we are now able to write 
down their corresponding RG equations 



u4-;T b R (p h q\ { ) =-u-^AT b R ( P[l , qil ) 



did 



x [lip\\^) +l(p\\ ~Q\ 
u-^T c R (p hqi] ) = -cj—AT R (p {l , qil ) - 



T b R (Pb<l\\) 
)] (4.11a) 



dio 
x [7CP|h u ' 



2^0||>«||) 



)+l{l\\-P\\M] (4-llb) 



where b = CDW, SDW and c = SSC, TSC. The full 
expressions for the symmetrized counterterms A7^ fl 's for 
the anomalous dimension is given in Appendix A. Due 
to the particular shape of our flat FS, the renormalized 
couplings must be symmetrical with respect to the ex- 
change of upper(right) and lower(left) particles and the 
change of sign of the external p,-ii 's. It then follows that 



■L CDW(SDW) CDW(SDW| 




R± 



sscfrsq 



R± 



ssctrsc) 



FIG. 7: The respective susceptibilities generated from the 



renormalized form factors T, 



R± 

CDW (SDW) 



and T* L± 



SSC(TSC)- 



rrR± 
2 b 



(PH , 9|| ) = T«(p„ , n )± T«(-p„ , gn ) (4.13a) 
T c fl± (p||,g||) =T R (p ll , q[l )±T R (-p ll , qil ) (4.13b) 



where again b = CDW, SDW and c = SSC, TSC. The 
(+) sign is associated with the s-wave symmetry whereas 
the (— ) sign is associated with the d-wave symmetry in- 
stead. Now, considering the symmetries of the couplings 
(|4.12|) and the new symmetrized form factors (|4.13a|) 
and (|4.13b|) we can write the RG equations (|4.11a|l and 
(|4.11b|) in the following way 



^ ± b. ll ,, ll ) = -^Ar-(p ll , 9|| ) + i^(p ll , 9|l) 

x [l(Pb^) +l(P\\ ~ Q\\,u)] (4.14a) 

^ ± b. ll , 9|| ) = -^Ar-(p ll , 9||) + i7;-(p ll , 9|l ) 

x biPb") +7(5|| -P\\,u)] (4.14b) 



R± 

CDW' 



-R± 



AT K± AT 1 



where the expressions for AT C 

/S.T R ^ C can be found in the Appendix B. The plus sign in 
the DW's is associated with the charge and spin density 
waves. However, the minus sign symmetry of the paral- 
lel momentum along the FS in the DW's yield circular 
charge(spin) currents flowing around the square lattice 
with alternating directions. In this way we associate this 
symmetry of the DW's with the charge and spin cur- 
rent waves also known as flux phases. Once the renor- 
malized TcDW(SDW) an d ^"sSC(TSC) are f° un( l we can 
define the related susceptibilities following the diagram- 
matic scheme shown in Fig. As one can see there is a 
IR divergent bubble in each channel. As a result we get 



X^fell,") 



4ir 2 v F ln (fi) 



T>3 



dp,, {T R± (jp hn )) 



*-r b R± (Pb<i\\ 



(4.15a) 



A-K 2 Vp 



xT c R± ( Pll , q{l ) (4.15b) 
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where T>-$ and T>± are the intervals determined in the Ap- 
pendix B and b and c refer to the symmetries mentioned 
before. Notice that if we don't go beyond one-loop or- 
der all RG equations for the susceptibilities approach the 
strong coupling regime. As a result all susceptibilities di- 
verge in the IR limit at that order of perturbation theory. 
In two-loop order we show that only one susceptibility 
at a time remains divergent and all the others approach 
fixed plateaus values instead. Having said that we also 
call attention to the fact that in two-loops the leading 
susceptibility approaches the strong coupling regime in a 
much slower rate than the one found in one-loop order. 
This might mean that the IR divergence may disappear 
altogether in the presence of higher order corrections and 
as a result one should not expect true long-range order if 
all this is taken into account. We discuss this further later 
on. From Eqs. i|4.15a|) and l|4.15bjl we arrive immediately 
at the RG equations for the various susceptibilities 

T b fl± (p||,g||) (4.16a) 

r c fi± (p||,g||) (4.16b) 

These expressions are written in general form. How- 
ever, since we are just considering gm and giR couplings 

it turns out that T^f = (V^) . 

V. RESULTS 

In order to solve all those RG equations we have to 
resort to numerical methods since we want to estimate 
how the order parameters change as we vary the scale lo 
to take the physical system towards the FS. As we have 



1.0 




Step I 

FIG. 9: Susceptibilities X^ b (si\\ = 0; I) against step 1 for one- 
loop approach with g 1R = g 2R = 10 as initial conditions for 
the couplings. The SDW s-type diverges for 1=1.2. 



done in an earlier work we discretize the FS continuum 
replacing the interval —A ^ pn ^ A by a discrete set 
of 33 points. For convenience, we use u> = f2exp(— I), 
where = 2kp\ with I being our RG step. We choose 
£l/vpA — 1. In view of our choice of points for the FS, 
we are only allowed to go up to I 2.8 in the RG flow to 
avoid the distance u> to the FS being smaller than the dis- 
tance between neighboring points since the discretization 
procedure no longer applies in this case. 

The RG equations for all couplings considered here 
were reproduced in section III. As a result all RG equa- 
tions (|XT|) . J33JI, H4.11a[) and (|4.11b|) and (|4.16a() and 
(|4.16b|) have to be solved simultaneously. Hence, the nu- 
merical estimates become much more involved. 

To improve the understanding of our results we also 

0.4t 1 
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consider the one-dimensional limits of all our RG equa- 
tions. We reproduce all those one-dimensional RG equa- 
tions in our appendix C. Following the same scheme used 
for D = 2 all one-dimensional RG equations are solved 
numerically in a self-consistent manner. Notice that as 
expected the equation ud(giR — 2g2R)/duj = is triv- 
ially satisfied. In Fig. [SJwe display the one-dimensional 
phase diagram in coupling space. This phase diagram is 
constructed in the usual way taking into account the two 
dominant susceptibilities that diverge for a given initial 
values of the renormalized couplings giR and g 2 R- We do 
that because there are always two divergent susceptibil- 
ities in all region of the phase space. Our numerical es- 
timates are in good agreement with Solyom's earlier RG 
results 13 . However, we noticed that those divergences 
take place only when the RG step is too large or, equiva- 
lently, u> is very small. This effect is amplified as we move 
from one-loop to two-loops. As a result despite the fact 
that there might be strong collective fluctuations in such 
a system we cannot say that they are of long range type 
since when we consider higher orders in our perturbation 
theory the divergences only take place in even larger RG 
steps. This situation will reflect itself more emphatically 
in the two-dimensional case as we will see next. 

Now, we move on to our original RG equations in two- 
dimensions. The choice of the initial conditions at I = 
in the RG equations are, in principle, arbitrary. However, 
as we have done before, we set the couplings equal to each 
other to mimic the Hubbard on-site repulsive interaction 
parameter U in our RG scheme. Namely g 1R — g 2 R = 10, 
where g iR = giR/itvp. In another paper we give more 
details about this connection with the Hubbard modeli^. 

In order to reproduce some symmetries of the order 
parameters with respect to FS we choose the following 



> 15 



initial conditions (I = 0) 



~ 12- 



Q. 

CD 
U 
» 
3 

w 

aT 
d. 
>. 

-4— • 

w 



Q 




1.0 1.5 2.0 
Step I 



r CDw(Pb1\\) 

T CBW(P\\:1\\) 



= t sdw (P|| 1 9|| ) = T ssc (P|| . 9|| ) = 

T TSc(Pb<l\\) = 1 



wjw(P||i«I|) = ^o(P||i9||) 



(5.1a) 



T TSc(PbQ\\) = ^ sin (ct) ( 5 ' lb ) 



Notice that despite the fact that the initial values of 
the form factor are either unity or a function of j> | the q\\ 
dependence is generated naturally by the renormaliztion 
process. The choice made for the symmetrized (+) ver- 
tices is motivated by their independence with respect to 
the change of sign of pii . In contrast with the antisym- 
metric vertices our choice is oriented by our need to re- 
produce the symmetries of the flux phases and the d x 2_ y 2 
superconductivity with respect to the FS. Furthermore, 
we take all susceptibilities equal to zero at I = 0, that is, 
x ± ( q]l ; / = 0) = (i = CDW, SDW, SSC, TSC) . 

To arrive at the final results we do a simultaneous cal- 
culation of the renormalized form factors and the corre- 
sponding susceptibilities together with the flow equations 
for the coupling functions g\R and g 2 R- It emerges from 
our numerical estimates that the SDW of symmetry type 
s produces the dominant susceptibility when the transfer 
momentum is equal to the nesting vector (q = Q*). To 
make the comparison of our results with previous esti- 
mates found in the literature we display our computa- 
tions in two steps. Initially we show the one-loop order 
results since they can be compared directly with calcu- 
lations presented by other groups. Next, we move on 
to display our two-loop calculations. Following that we 
discuss the most interesting results and the difference be- 
tween our estimates in one-loop and two-loops. 




FIG. 11: SDW+ susceptibility Xsdw(i\\ > with g 1R = g 2R = 
10 as initial conditions for the couplings. The SDW s-type 
diverges for I = 2.15. Dashed line refers to perfect nesting of 
the FS. Solid line refers to the RG step I — l c utoff considering 
a "corrugation" in the FS. 



FIG. 12: SDW+ susceptibility X £ 

Dw(l\\' against Q\\ for sev- 
eral RG steps with g 1R — ~g 2 R = 10 as initial conditions for 
the couplings. 
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FIG. 13: The leading susceptibilities for several initial values 
of coupling functions. 

A. One-loop RG approach 

The one-loop results for all susceptibilities are dis- 
played in Fig. |5| As can be seen all susceptibilities di- 
verge but the leading one is the SDW+ symmetry in the 
repulsive HM like regime. According to this, we should 
expect an insulating spin density wave state and no sign 
of even nonconventional metallic behavior in the physical 
system. It is interesting to observe that the second most 
pronounced renormalized susceptibility corresponds to a 
d x 2_ y 2 superconductivity symmetry. 

Our result is consistent with other approaches based 
on one-loop fermionic functional RG - and the so-called 
parquet methodii. In all those approaches, the feedback 
of the quasiparticle weight Z is not taken into account. 
The fact that Z — > produces drastic changes in the flow 
of all RG equations in two-loops as we will see next. 

B. The full two-loop approach 

In this final case, we solve the Eqs. I|4.11a|) and (|4.11b|) 
and Eqs. I|4.16a|l and (|4.16b|) simultaneously with the RG 
equations for the renormalized couplings using the same 
initial conditions given by Eqs. (|5.1a|) and l|5.1b|) . namely, 
xt(<l\\ 'J = 0) = and g 1R = g 2R = 10 respectively. 

The results obtained for all susceptibilities excluding 
the Xsdw+ are shown in Fig. ^| As one can see the 
effect of Z in the RG equations changes drastically the 
one-loop scenario. Notice that the plateaus values ap- 
pear at I = 2.2 for almost all susceptibilities, that is, 
when the strong renormalization of the coupling func- 
tions takes place. In contrast with the picture described 
for the coupling functions at the same two-loop order, 
these plateaus are real fixed points since their values are 
not sensitive to our FS discretization procedure. Thus, 
we can infer from this, that with the exception of the 
SDW+ symmetry, all existing order parameters appear 



as short range correlations. However, for this HM like 
scenario, the SDW+ symmetry manifests itself in a com- 
pletely different manner and to illustrate this we plotted 
the SDW+ susceptibility as a function of the parallel 
transfer momentum (q») in Fig. 1121 

As g|| approaches zero the SDW+ susceptibility in- 
creases drastically as a function of the RG step I. In 
contrast for any other value for the parallel transfer mo- 
mentum (<7|| 7^ 0) the RG equation for Xsdw fl° ws to a 
given fixed value. Does this means that we can assign 
to the SDW+ divergence at qu = a long-range type 
correlation? As we emphasized earlier the effect of the 
quasiparticle weight in two-loop approach the SDW+ 
divergence becomes predominant for larger values of the 
renormalization step I. This is a strong evidence that 
even this divergence might not be what it seems. In 
fact, we expect that if we were to consider higher or- 
der contributions in perturbation theory this divergence 
would occur for even larger I values. In order to interpret 
correctly the meaning of this divergent susceptibility we 
have to keep in mind the limitations of our RG analy- 
sis. From the start we said that we were dealing with a 
two-dimensional flat FS with rounded corners. This FS 
represents a doped electron system slightly away from 
half-filling. As a result there exists a non-zero chemical 
potential \x which is an indirect measure of that doping 
regime, [i introduces a new scale into the problem. At 
RG scales ui » fj, we can safely treat the FS as being 
made of perfectly flat patches with rounded corners. As 
lu is lowered (i.e. as we increase the number of RG 1 
steps) our RG resolution of FS must be corrected to in- 
clude its "corrugation" effects. To take this resolution 
effect approximately into account we follow Zheleznyak 
et aiii in establishing a cutoff l c beyond which we are no 
longer able to be sure about the perfect flatness of our 
FS. 

We take /x such that l c = ln(Sl cuto ///|/i|) = 1.95 since 
Xsdw OIu y diverges for I — 2.15. This means that for 
I = l c the SDW s-type will now reach a plateau value 
before this divergence takes place. We display this new 
Xsdw m Figure ITT1 This excludes the existence of long- 
range order in such a lightly doped Hubbard model like 
regime. 

In Figure IT^l we display a phase diagram with the lead- 
ing susceptibilities for several initial conditions of the 
coupling functions. In our calculations we notice that 
when we set initial values for the couplings nearby each 
axis the symmetry associated with the next phase be- 
comes more pronounced. As an example, if we are in the 
SDW+ phase and we choose initial conditions nearby 
the straight line g 1 = 2g 2 we obtain as a result a increase 
in the TSC+ susceptibility and so on. Moreover, we 
can emphasize that this phase diagram is quite different 
from that in one-dimension. The latter shows the pos- 
sible ground states inferred from the renormalized form 
factors which in turn are a mixture of two states while in 
two-dimensions there is only one pronounced susceptibil- 
ity at a time which we associate with a strong collective 
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oscillation. This shows how much the one-dimension sce- 
nario can change as one goes from one to two-dimensions. 



VI. CONCLUSION 

In this paper we have performed a full two- 
loop theoretical RG calculation for the symmetries 
CDW±,SDW±,SSC± and TSC± in the presence of 
a square two dimensional FS with rounded corners. We 
neglected Umklapp effects at this stage since we are not 
at half-filling. To render the theory finite we applied the 
RG field theory by adding appropriate counterterms or- 
der by order in perturbation theory. 

Due to the particular shape of our FS we can re- 
produce some symmetries by just choosing appropri- 
ate initial conditions for the renormalized form factors. 
The general expectation is that the resulting physical 
state has several competing instabilities such as the flux 
phases (CDW— and SDW— ), the singlet d x 2_ y 2 super- 
conductivity (SSC-), the d-type triplet superconductiv- 
ity (TSC-), the singlet and triplet superconductivity of 
s-type (SSC+ and TSC+) and finally the spin density 
wave (SDW+) and charge density wave (CDW+). At 
one-loop order all those symmetries diverge and the lead- 
ing divergence is SDW+ suggesting that there is a long- 
range type correlation. 

In two-loop order for the repulsive HM like initial con- 
ditions we obtained that all susceptibilities go to stable 
fixed points except the SDW+ symmetry that continues 
to diverge in this case. In order to interpret correctly the 
meaning of this divergent susceptibility we call attention 
to the fact that the divergence in two-loop order takes 
place for larger values of the RG step I in comparison 
with what is obtained in one-loop order. It is therefore 
reasonable to expect that in higher orders the divergence 
takes place for even higher values of I, This happens due 
to the limitation of the RG analysis. 

In this model we assumed that the FS is completely 
flat and, as a result, perfectly nested. However, in real 
situations, this naive assumption breaks down since most 
FS's are expected to exhibit a finite curvature at some 
critical scale uj c as one approaches the low-energy limit. 
In such a condition the nesting is destroyed. If this hap- 
pens, the RG flow of the susceptibility associated with the 
SDW+ order parameter will become interrupted due to 
the fact that this quantity strongly relies on the nesting 
properties of our FS model. Consequently, this result as- 
serts that there should be only short-range order in the 
system in agreement with the general expectation that 
strong quantum fluctuations should be a dominant fea- 
ture in such low dimensional systems. 

We sketched the phase diagram for different renor- 
malized coupling constants initial conditions in two- 
dimensions. There is only one pronounced susceptibility 
at a time which we associate with a dominant collective 
oscillation. In order to test our method in well known 
grounds we solved all RG equations for one-dimensional 



case. We reproduce correctly the phase diagram dis- 
played earlier by Solyomi^ together with the possible 
ground states of the system inferred by renormalized form 
factors. 

To have a better idea about the nature of the domi- 
nant instability we have to calculate the corresponding 
uniform susceptibilities which are calculated taking the 
transfer momentum equal to zero (q = (0,0)). This has 
been done in a recent paper for the particular case of 
the repulsive HM like initial conditions . We obtained in 
case that the charge compressibility flows to zero, that is, 
dn/dfi — > and, hence, we could indeed consider that the 
SDW+ divergence would lead to a symmetry breaking 
contradicting all arguments presented previously. How- 
ever, another quantity flows to zero as well. When we cal- 
culate the uniform spin susceptibility we find the opening 

r ■ uniform n < TT • 

ot a spin gap as X S pin flows to zero. Hence, it is our 
belief that the quantum regime associated with such a 
state is characteristic of an insulating spin liquid. This 
conclusion is in agreement with the results presented in 
this work. 

This work was partially supported by the Conselho 
Nacional de Desenvolvimento Cientfico e Tecnolgico 
(CNPq). 



APPENDIX A 

In this appendix, we write down the explicit form of the 
ATfs (i = CDW, SDW, SSC, TSC) a n d 7 's w h ich are 
taken into account in the Eqs. I|4.11ajl . (|4. 1 1 b|) . (|4.14a|) 
and (|4.14b|) . We also give the several intervals of inte- 
gration that are considered throughout this work. They 
are following 




-A k\\ ^ A, 
-A p|| s$ A, 
-2A < <7|| sC 2A, 
-A p|| - g|| A. 



{-A k\\ A, 
-A p|| s$ A, 
-2A sC g|| 2A, 
-A < g,| - fc|| < A, 
-A < 9|| -p|| < A. 



r a < pi, < a, 

V 3 = I -2A s? 9|| < 2A, 

[ -A<j»|| -9|| < A. 



r a < p|| < a, 

P 4 = I -2A sC 9|| 2A, 
[ -A < 9|| -p|| < A. 
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2?B = 



V* = 



-A jfe|| A, 

-A<p|| s; A, 

-A ^ ?! || < A, 

-A -fc|| +p|| + q n A. 



-A < fey < A, 
-A < p|| s: A, 
-2A < 3|| < 2A, 
-A < q n < A, 

-A < -k\\ + p\\ - q\\ + qi\\ < A. 



k\\, )9iR {p\\, +2g2 R {p\\,qi\\,-k\\ 
x-92R(k\\,-k\\ +p|| -.gifl(p||,9i||,fc||) 

X 92R ( fc |h~ fc || + P|| + Si 11,3111) ~ 92R{p\\-,<li\\, 

-fc|| +gi||)pifl + 3i||,P||)L (A5) 

APPENDIX B 



V- 



-A< *n < A, 
-2A < g|| 2A, 
-A<p|| < A, 
-A<g 1( | < A, 
-A -fc|| +?|| -p|| 



- S11 



< A. 



We begin with the expressions for the AT^'s associ- 
ated with r( 2 '^'s. We get 



AT/ 1 



cdw 



(Pb1\\) = 



rW 2 



x / [2g 1R (k\\,p\\ - <7||,P||) 



,P\\ -Q\\,h ~ 1\\)] T CDw( k h<l\\), 



g2R 1 



(Al) 



Considering the symmetries obeyed by the coupling 
functions we get the symmetrized (±) renormalized form 
factors whose counterterms A7J are given by 



AT R± 



(P||>9||) 



x / dfc|| [2g 1R {k\\,p\\ - q\\,P\\) 

~g2R (k h P\\ -Qh k \\ ~ 1\\)] r CDw{ k \\^\\)^ 



AT R± 



tW 2 



(Bl) 



x / dk ] \g 2R (k h p l \ -q\\,k\\ - q\\) T^ DW (k\\,q\\), (B2) 

J-Dt 



AT? 



SDW 



(Ph1\\) 



x / dk\\g2R{k\\,pi\-qi\,k\\-q\\)T SDW (k\\,q\\), (A2) 



AT ssc(Pb<l\\) 



x / dk\\[g 1R (k\\,-k\\ +<?||,P||) 

+ g2R (k\\,-k\\ +q\\,-p\\ +q\\)]Tssc(h>q\\)> ( A3 ) 

and finally 



AT R± 



(Plh^ll) 



4:Tt z vf \u 



x / dk\\[gi R (k\\,-k\i +q\\,p\\) 
Jt> 2 

+ g 2R (k l] ,-k ]l +q l] ,-p ll + q ll )]T*±(k ll ,q ]l ), (B3) 



and finally 



x / dk\\ [gm (k\\,-k\\ + <7||,P||) 
Jt> 2 

-92R {k\\i~ k \\ + 1b~P\\ + ( l\\)] T rsb( k \\^\\)- 



(B4) 



AT R 



TSC 



(P||.9||) 



x / dfty [.gifi (fc||, + q\\,P\\) 



g 2 R{k\\,-k\\+q\\,-p\ ] +q{)]T^ sc {k h q\ ] ). (A4) 



APPENDIX C 

In this appendix we will present all RG equations up 
to two-loop for the one-dimension case. For the quasi- 
particle weight we get simply 



The anomalous dimension 7 used in our RG equations 
for the renormalized form factors is given byi^ 



UJ- 



d\nZ 

dui 



4ir 2 v F 



(9ir + 92R ~ 9iRg2R) = 7 (CI) 
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In ID the flow equations for the renormalized coupling 
functions reduce simply to 



dgw _ g( R . g{ R 



duj ttvf 2ttv 2 



(C2) 



F 



dTsscR 1 



w- 



2 3 duj 2ttvf 



-(giR + 92r)Tsscr + iTsscr (C6) 



dg2R _ gl R , gf R 



dT^rscR 1 

( C3 ) ^ — do; — 1 = 2ttv^^ 92R ~ gw ^ TSCR + "(^tscr (C7) 



Finally, the renormalized form factors are much simpler 
and they are now given by 



Consequently, the ID resulting susceptibilities are now 
determined by 



dTcDWR 1 , ,rr , rr 

V : — ~ K^giR — g2R.)J-CDWR + l-^CDWR 



duj 2ttvf 



(C4) 



dXal 



duj 2-kvf 



' (%r) 2 



(C8) 



dTsDWR -g2R rr fnK , 
W j = 7. -iSDWR + 7-lSDWR (^O) 

dui Zttv f 



where a = CDW, SDW, SSC, TSC. 
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